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ABSTRACT: A least squares regression interval mapping model was derived to detect quantitative trait loci (QTL) with a unique 
mode of genomic imprinting, polar overdominance (POD), under a breed cross design model in outbred manmials. Tests to differentiate 
POD QTL from Mendelian, paternal or maternal expression QTL were also developed. To evaluate the power of the POD models and to 
determine the ability to differentiate POD from non-POD QTL, phenotypic data, marker data and a biallelic QTL were simulated on 512 
F2 offspring. When tests for Mendelian versus parent-of-origin expression were performed, most POD QTL were classified as partially 
imprinted QTL. The application of the series of POD tests showed that more than 90% and 80% of medium and small POD QTL were 
declared as POD type. However, when breed-origin alleles were segregating in the grand parental breeds, the proportion of declared 
POD QTL decreased, which was more pronounced in a mating design with a small number of parents (Fq and Fi). Non-POD QTL, i.e. 
with Mendelian or parent-of-origin expression (complete imprinting) inheritance, were well classified (>90%) as non-POD QTL, except 
for QTL with small effects and paternal or maternal expression in the design with a small number of parents, for which spurious POD 
QTL were declared. (Key Words: QTL, Imprinting, Polar Overdominance, Pig, Breed-cross) 



INTRODUCTION 

Genomic imprinting, a non-Mendelian form of gene 
expression, has been one of the interesting research areas in 
mammahan genetics due to its unique expression and 
inheritance mechanisms and its important role in growth, 
development, and behavior (Tycko and Morison, 2002; 
Wilkens and Haig, 2003). One typical imprinted gene, 
insulin-like growth factor 2 {IGF2), has been shown to 
influence body composition in swine (Jeon et al., 1999; 
Nezer et al., 1999). Genomic regions with imprinted effects 
can be detected using QTL interval mapping in F2 crosses 
between lines or breeds that segregate for marker alleles 
because of the abihty to follow parental origin (Thomsen et 
al., 2004; Choi et al., 2011; Li et al., 2011). Development 
and appUcation of statistical methods to detect 
chromosomal regions or quantitative trait loci (QTL) with 
imprinting effects in a breed cross design have been 
developed and implemented for F2 crosses of breeds of 
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swine (Knott et al., 1998; de Koning et al., 2000; Ruckert et 
al., 2012). A specific mode of genomic imprinting that has 
been uncovered in mammals is polar overdominance (POD), 
for which one type of heterozygous progeny differs in 
phenotype from the other heterozygous type, depending on 
the parent from which the mutant allele was received, and 
also differs from the two types of homozygotes. The only 
known phenotype caused by this unique expression in 
mammals is the 'callipyge' phenotype, which is expressed 
as muscular hypertrophy in sheep (Cockett et al., 1996). 
Subsequent evaluations found that the lambs with the 
callipyge-causing genotype at the callipyge locus (CLPG) 
had heavier carcass, greater fat-free lean, and less fat than 
those with the other genotypes (Freking et al., 1998). In 
swine, a DLK polymorphism in the CLPG gene was also 
reported to be associated with growth and fatness in a cross 
of commercial pig breeds (Kim et al., 2004). 

The first objective of this study was to derive a POD 
model for a genome-scan in a cross between outbred breeds 
or lines, so as to enable a systematic search of other regions 
or loci that exhibit POD across the genomes of outbred 
animals. The second was to vaUdate the proposed POD 
models by performing simulation studies with different 
mapping designs and sample sizes, and different modes of 
inheritance, QTL effects, and allele frequencies in parental 
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breeds. 

MATERIALS AND METHODS 

Polar overdominance QTL models 

Derivation of the model was based on the framework of 
the least-squares regression interval mapping method for 
detection of QTL in a breed-cross model, which detects 
differences in the average effects of alleles that originated 
from the two parental breeds (Haley et al., 1994). There are 
four possible QTL genotypes among F2 animals in the 
breed-cross model based on breed origin, 11, 12, 21, 22, (1 
and 2 represent alternate QTL alleles derived from parental 
breeds 1 and 2, and the first allele represents the allele 
received from the Fi sire). Following de Koning et al., 
(2002), the full imprinting model for a putative QTL at a 
given position is: 



animal at the putative QTL position conditional on flanking 
marker genotypes, such that Ppat = (Pii+Pi2)-(P2i+P22) and 
Pmat = (Pii+P2i)-(Pi2+P22), whcrc Py is the vector of 
probabilities of inheriting QTL genotype ij. The effects of 
the four possible QTL genotypes in terms of parameters of 
the full model (1) are (apat+a^at), (apa-amat+d), (-apat+a^at+d), 
and (-apat-a,nat), respectively (Figure 1). To generate POD 
models, restrictions were put on the QTL genotype effects, 
such that in a type I POD model (here we refer to type I and 
type II POD when the differential phenotype is observed for 
QTL genotype 12 and for genotype 21, respectively), 
genotype effects of 11, 21, and 22 are the same, such that 
(apat+amat) = (-apat+amat+d) and (-apat+a^iat+d) = (-apat-a^iat). 
This makes apat = l/2d and a^jat = -l/2d, where d represents 
the dominance effect of the 12 genotype over the average of 
the three other genotypes. Using this parameterization in 
model 1, the type I POD model becomes: 



Y - Xb-napatPpat+amatPmat+dPd+e 



(1) 



Where, Y is a vector of phenotypes on F2 animals, X is 
a design matrix, b is a vector of fixed and covariate effects, 
and e is a vector of uncorrected residuals with constant 
variance. The coefficients apat and a,nat are the effects 
associated with the paternally and maternally inherited 
breed specific QTL alleles, respectively, and Ppat and P^jat 
are the corresponding vectors with probabilities for each 



Y = Xb+d(l/2Ppa,-l/2P„a,+Pd)+e 
which reduces to: 

Y = XbH-2dPi2+e 

Similarly, the type II POD model is: 

Y = XbH-2dP2i+e 
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Figure 1. Models and coefficients for QTL effects for each genotype in a cross between two outbred breeds or lines with alternate QTL 
alleles (1 and 2) in a three generation breed cross population. In the full model, both homozygotes have different additive effects 
depending on inheritance of alternate alleles through the Fi sire (apat) and Fi dam (a,^a,). Heterozygotes have an additional dominance 
effect (d). Restrictions (in parentheses) were imposed on QTL effects in the full model (De Koning et al., 2002) to identify: Mendelian 
QTL when there is no difference of additive effects between alleles derived through the F] sire and Fj dam; paternally (maternally) 
expressed QTL, when alternate alleles inherited through Fi sire (dam) differ in effect; type I (II) polar overdominance QTL, when the 
differential phenotype is observed for QTL genotype 12 (21). 
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Parameterizations to derive Mendelian, paternal, and 
maternal expression models are in Figure 1, following de 
Koning et al. (2002), resulting in: 

Y = Xb+aPa+dPd+e (4) 

Y = Xb+apatPpat+e (5) 

Y = Xb+a^,Pmat+e (6) 

Where a is the additive QTL effect and Pa is a vector of 
probabilities (P11-P22), and the other variables are as in (1). 

To estimate genetic variance due to a gene with POD 
effect, the population mean at the QTL under random 
mating was derived with parameters p and q (alternate allele 
frequencies from the two parental breeds), and d 
(dominance effect as in the above POD model). If a gene 
has complete POD inheritance mode, then one heterozygote 
has a genotype effect of 2d and the other three genotypes 
have no additive (a) or dominance effect. The population 
mean then becomes 2dpq and the genetic variance due to 
the POD QTL is 4pqd\l-pq). Under random mating, 
breed-origin allele frequencies, p and q are equal to 0.5 in 
the F2S of a breed-cross population, such that the mean and 
variance become 0.5d and Q.lSd^. For Mendelian and 
complete imprinting QTL, genetic variances are 2pq[a+ 
diq-p)f+i2pqdf and 4pqa^, respectively (de Koning et al., 
2002). In an F2 population where breed-origin QTL alleles 
have equal frequency (p = q = 0.5), the variances become 
0.5a^-H0.25<i^ and a^, respectively. 

To detect POD QTL, an F-test statistic for models 2 and 
3 against the null hypothesis of no QTL (Fpoo) can be 
obtained at unit increments along the chromosome (Test 1). 
For significant QTL based on Test 1, the following two sets 
of tests were performed to determine whether the POD 
model fits better than the full, Mendelian, paternal, or 
maternal expression models (Tests 2 and 3). 

Test 2. Flof- Ha = full imprinting model (1) vs Ho = 
POD model (2 or 3). This test was conducted at the most 
likely position under the POD model to determine whether 
the POD model does not have a significantly poorer fit than 
the full imprinting model. Test statistic F^of has an F- 
distribution with 2 and n-4 (n is the number of F2S) degrees 
of freedom. 

Test 3. Ha = POD model (2 or 3) vs Ho = MendeUan (4), 
paternal (5), and maternal (6) expression models, using the 
following likelihood ratio test statistics, respectively: 

LRTmend = n * ln(RSSMEisnyRSSpoD)> 

LRTpAT = n* ln(RSSpAT/RSSpoD), 



LRTmat = n * ln(RSSMAT/RSSpoD), 

Where, RSS is residual sum of squares of the 
corresponding model at the most likely position under the 
POD model, and n is the number of individuals in the F2. 
Since these are not tests between nested models, 
significance thresholds were derived empirically by 
randomly switching elements of P vectors (within 
individual); Pn vs P12 for LRTpat and P12 vs P22 for LRTmat 
in the type I POD model; P21 vs P22 for LRTpAx and Pn vs 
P21 for LRTmat in the type II POD model; and P12 vs P21 for 
LRTmend in both models. A QTL with POD mode of gene 
action was declared when Ho wa not rejected in Test 2 and 
all Ho's were rejected in Test 3. 

Mendelian vs parent-of-origin expression analysis 

POD QTL were tested with a series of tests to 
differentiate parent-of-origin effects (complete imprinting) 
from Mendelian effects. F-tests between Mendelian and 
parent-of-origin models were applied according to the 
decision tree described in Thomsen et al. (2004), with some 
minor modifications of the specific tests conducted. First, 
the Mendelian model (4) was tested against the null model. 
If significant at a 5% chromosome- wise (CW) level, the full 
model (1) was tested against the Mendelian model at the 
best position of the full model within the QTL region. If not 
significant at the 5% CW level, a Mendelian QTL was 
declared, otherwise a QTL with parent-of-origin effects. 
When the Mendelian model was not significant in the first 
step, the full model was tested against the null model and a 
QTL with parent-of-origin effects was declared if 
significant at the 5% CW level. To determine the nature of 
the detected parent-of-origin effects, tests of the full against 
the paternal (5) and maternal (6) models were conducted at 
the best positions of these respective models within the 
QTL region. If the full model explained significantly more 
variation than the paternal model (at the 5% CW level) and 
the full model did not explain more variation than the 
maternal model (5% CW), then fl(pat) = 0 and a maternally 
expressed QTL was inferred. A paternally expressed QTL 
was inferred vice versa (a(mat) = 0). If tests of the full against 
the paternal and maternal expression models were both 
significant or both not significant, partial expression was 
inferred (Thomsen et al., 2004). Chromosome-wise 
thresholds for the tests between alternate expression models 
were set equal to thresholds obtained for tests against the 
null hypothesis of no QTL with equivalent degrees of 
freedom in the numerator of F-test statistic, as suggested by 
Thomsen et al. (2004). Thus, thresholds for tests of the 
paternal against the null model were used for tests of the 
full against the Mendelian model, and thresholds for tests of 
the Mendelian against the null model were used to test the 
full model against the paternal and maternal models. 
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Simulation 

To validate the POD models and tests and to evaluate 
their power and ability to distinguish alternative QTL types, 
F2 populations were simulated based on two designs that are 
relevant to experimental pig QTL mapping populations 
(Bidanel and Rothschild, 2002). Design I mimicked the 
Berkshire and Yorkshire F2 cross described in Malek et al. 
(2001). This design was based on a Umited number of 
parents, with two Fq grand sires of one breed and ten Fq 
grand dams of another breed, to generate 10 Fi's per dam. 
Eight Fi sires and 32 Fi dams were randomly chosen to 
produce 16 F2's per Fi dam, for a total of 512 F2 offspring. 
The second design was based on a larger number of parents, 
with 20 Fq sires and 80 Fq dams, with five offspring per dam. 
A total of 19 Fi sires and 57 Fi dams were randomly chosen 
to produce 513 F2 progeny (9 per full- sib family). A 
chromosome of 100 cM was simulated with 11 markers at 
10 cM intervals. Markers were simulated with four alleles 
with frequencies of 0.6, 0.2, 0.1, 0.1 in one breed and 0.1, 
0.1, 0.2, 0.6 in the other breed, following simulations that 
were designed to mimic the marker information content of 
the Berkshire and Yorkshire cross (Malek et al., 2001). 
Biallelic QTL with different modes of inheritance pattern 
were simulated at position 75 cM on the chromosome: 
additive, complete dominance, paternal and maternal 
expression, and type I and II POD effects. QTL effects were 
defined as large, medium or small QTL, such that the QTL 
explained 32%, 12.5%, or 5.1%, respectively, of the 
phenotypic variance for the quantitative trait that was 
simulated. Different QTL effects were set depending on 
inheritance pattern and defined QTL size, such that the 
defined QTL sizes were consistent across QTL types (Table 
1). Alternate QTL alleles in parental breeds were assumed 
to be fixed (1.0/0.0) or to be distributed with different 
frequencies (0.8/0.2). Five hundred replicate data sets were 
simulated for each set. Thresholds at the 5% CW level for 
QTL detection in each QTL model were derived from three 
thousand replicates with QTL effects set to zero (Test 1). 
Significance thresholds for Tests 2 and 3 were determined at 
a comparison-wise 5% level, p Values for Test 2 were 
obtained using standard F statistic thresholds. Significance 
threshold values for LRTs (Test 3) to test QTL with POD 
inheritance versus Mendelian, paternal, and maternal 



inheritance were obtained with 2000 permutations for each 
replicate. 

RESULTS 

Detection of POD QTL 

Table 2 presents power to detect POD QTL under 
different genetic models, estimates of QTL position and 
effects, and of the proportion of phenotypic variance 
explained by POD QTL under different mapping designs, 
QTL types and sizes, and allele frequencies in parental 
breeds. Most POD QTL were detected in their 
corresponding POD models when alternate alleles were 
fixed in parental breeds, for all QTL sizes, i.e. power to 
detect QTL was close to 100%. Also, POD QTL with large 
or medium size was detected in imprinting models at 
similar detection rates as in the POD models. However, 
when QTL alleles were segregating in the parental breeds, 
power to detect a POD QTL using their respective POD 
models decreased, especially for small QTL in design I, 
where a small number of parental animals were used. When 
increasing the number of F2 animals (double the F2 size), 
more POD QTL were detected using their respective POD 
models when the allele segregated in the parental breeds, 
e.g. more than 10% increase for small QTL. 

The QTL that were detected with statistical significance 
under POD models were further classified as POD QTL, 
using the tests (2 and 3) described above. Generally, POD 
QTL were most often declared to be POD QTL when 
parental breeds were fixed for alternate alleles, e.g. more 
than 95% for large QTL, an around 90% and 80% for 
medium and small QTL, respectively. When QTL size was 
small, around 10% of POD QTL were declared to be type I 
or n (Table 2). When alleles were segregating in parental 
breeds, the proportion of declared POD QTL decreased 
substantially across POD QTL types and sizes, and mating 
designs. The degree of failure to declare a POD QTL was 
larger for design I than for design II, for large or medium 
QTL. When increasing the number of F2 animals, a 
marginal increase in frequency of declared POD QTL was 
observed, especially for the small QTL in design II (Table 
2). 



Table 1. QTL effect and variance for different QTL sizes and inheritance patterns in the simulation study 

QTL effect" 

QTL inheritance mode Genetic variance — 



Large Medium Small 



Additive (a) 0.5a^ 0.800 0.500 0.320 

Complete dominance (a = 0.5a^+0.25d^ 0.654 0.408 0.261 

Parental-origin expression (a) 0.566 0.354 0.226 

Polar overdominance (rf) OJScf 0.654 0.408 0.261 

" Different QTL effects under varying genetic models were set such that large, medium and small QTL explained 32%, 12.5% or 5.1%, respectively, of the 

phenotypic variance. Error variances were set 0.680, 0.875 and 0.949 for large, medium and small QTL, respectively, such that overall phenotypic 

variances become standard unit (1.0). 
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Table 2. Power to detect polar overdominance (POD) QTL under different genetic models, estimates and standard deviations (in bracket) 
of POD QTL effects and positions, and proportion of phenot5'pic variance explained by POD QTL based on simulated data 
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9 


30 


27 


23 


55 


8 


37 


6 


5 


6 


5 


29 


72.0 (16.8) 


0.19 (0.05) 


2.9 


(1.1) 


Design II (N = 


1026) 
































10.41 (1.0) 


100 


100 


100 


100 


97 


94 


0 


2 


0 


0 


0 


100 


74.9 (1.8) 


0.41 (0.04) 


12.6 


(2.2) 


10.41 (0.8) 


41 


84 


82 


96 


81 


65 


1 


9 


1 


13 


12 


71 


74.8 (5.6) 


0.22 (0.05) 


4.0 


(1.9) 


10.26(1.0) 


79 


92 


91 


100 


69 


86 


0 


8 


0 


7 


8 


85 


75.3 (3.8) 


0.26 (0.04) 


5.2 


(1.6) 


10.26(0.8) 


19 


53 


53 


84 


39 


60 


5 


7 


4 


14 


15 


42 


74.1 (9.8) 


0.16 (0.04) 


2.1 


(0.9) 


110.41 (1.0) 


100 


100 


100 


97 


100 


0 


96 


1 


0 


0 


0 


100 


75.1 (1.6) 


0.41 (0.04) 


12.6 


(2.3) 


110.41 (0.8) 


36 


80 


85 


78 


97 


3 


68 


10 


0 


12 


12 


72 


74.7 (6.9) 


0.22 (0.05) 


3.8 


(1.8) 


110.26(1.0) 


82 


94 


93 


67 


100 


0 


88 


6 


0 


6 


6 


89 


74.9 (4.2) 


0.26 (0.04) 


5.3 


(1.5) 


n 0.26 (0.8) 


17 


49 


51 


43 


84 


5 


58 


10 


4 


11 


11 


45 


73.8(13.1) 


0.16 (0.04) 


2.0 


(0.9) 



Polar overdominance (POD) QTL were simulated based on diiferent allele frequencies, magnitudes, and number of F2 progeny in two mating designs I: 2 Fq sires, 10 Fo dams, 8 
Fi sires for 512 or 1,024 F2 progeny; II: 20 Fq sires, 80 Fq dams, 19 Fi sires for 513 or 1,026 F2 progeny. A biallelic QTL was simulated at 75 cM for a 100 cM chromosome 
with 11 equidistant marlcers. Each marker had four alleles with different allele frequencies in parental breeds (0.6 (0.1), 0.2 (0.1), 0.1 (0.2), 0.1 (0.6) in breed A (B)). A total of 
500 repMcates were generated per each parameter set. 

* POD I or n refer to QTL, for which differential phenotype is observed for QTL genotype 12 (I) against genotype effects of U, 21, and 22, or genotype 21 (II) against 11, 12, 
and 22. POD QTL effects (d), 0.65, 0.41, and 0.26 were defined as large, medium or small QTL, such that the QTL explained 32%, 12.5% or 5.1%, respectively of the 
phenotypic variance. Alternate QTL alleles were homogenously (1 .0/0.0) or differently (0.8/0.2) distributed in Fq parental breeds. 

Proportion of replicates in which a POD QTL was detected at a 5% chromosome-wise level in Mendelian, paternal expression, maternal expression, POD I, or POD n models. 
If a POD QTL was detected at a 5 % chromosome-wise level in its respective POD model (Testl ), a series of POD tests (Tests 2 and 3) were performed for the POD QTL to be 
declared as POD I, POD II, or both types. 

'' POD QTL were tested with a series of tests to differentiate parent-of-origin effects (paternal, maternal, or partial expression) from Mendelian effects, according to THOMSEN 
etal. (2004). 

° Mean estimates (standard deviations) of the most likely QTL position from the replicates with POD QTL evidence in respective POD models. 

' Mean estimates (standard deviations) of POD effect (d) at the most likely QTL position from the replicates with POD QTL evidence in respective POD models. 

* Mean estimates (standard deviations) of proportion of phenotypic variance due to POD QTL (= 0.75d^) from the replicates with POD QTL evidence in respective POD 
models. 
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POD QTL were also tested for classification as non- 
POD inheritance types, i.e. MendeUan, paternal or maternal 
expression (Table 2). Most POD QTL were declared to be 
partially expressed for large or medium QTL when the 
frequency of alternate alleles in parental breeds was equal. 
When the QTL size was small, or alternate alleles were 
segregating in parental breeds, the proportion of POD QTL 
declared to be partially expressed QTL was substantially 
decreased. Instead, the proportion of POD QTL declared as 
paternally or maternally expressed QTL increased, up to 
around 35% for small QTL when alternate alleles were 
fixed in parental breeds. These POD QTL that were 
declared as parental origin expression QTL had a similar 
proportion of paternally and maternally expressed QTL 
across QTL sizes, POD QTL types (I, II) and mating 
designs (Table 2). With a larger number of F2S, the 
proportion of POD QTL declared as partially expressed 
increased, up to 32% (57% to 89% for small type II POD 
QTL in design II). Less than 10% of POD QTL were 
declared as Mendelian expressed QTL across all QTL sizes 
and types in both designs (Table 2). 

Estimates of POD QTL position and effect 

Position estimates for POD QTL that were detected in 
their POD models were close to unbiased and had high 
precision (low standard deviation) when power of QTL 
detection was high, i.e. for the large or medium QTL when 
alternate alleles were fixed in the parental breeds. When the 
QTL effect was small and alleles were segregating in the 
parental breeds, accuracy and precision of position 
estimates was the lowest (Table 2). Increasing sample size 
(number of F2S) gave more accurate and precise position 
estimates for a given QTL size and allele frequencies, 
which was consistent with power of detecting QTL. 

Estimates of POD QTL effects were also unbiased and 
precise for large and medium QTL with alleles fixed for 
alternate alleles in the parental breeds. However, when 
alleles were not fixed in the grand-parents, estimates of 
QTL effects were downward biased and had lower precision 
(larger standard deviation). Generally these estimates were 
obtained as a function of the absolute difference between 
the frequencies of alternate alleles in the Fq sire and dam 
breeds, i.e. 10.8-0.21 = 0.6, such that expected POD effects 
were 0.39, 0.25, and 0.16 for large, medium and small POD 
QTL, respectively. A larger number of F2 progeny gave 
more accurate and precise estimates of QTL effects, 
especially for small QTL with segregating alleles in the 
parental breeds (Table 2). Estimates of the proportion of 
phenotype variance explained by POD QTL were generally 
unbiased for large and medium QTL, slightly upward 
biased for small QTL, or downward biased for QTL with 
segregating alleles in the grand parents. 



Declaration of non-POD QTL as POD QTL 

Table 3 presents power to detect non-POD QTL under 
different genetic models, declared POD types in different 
mapping designs, QTL sizes, and different allele 
frequencies in the grand parental breeds. The power to 
detect QTL in the POD models was the greatest for paternal 
or maternal expression QTL (close to 100% for large and 
medium effect QTL with fixation of alternate alleles in the 
grand sire and dam parents), followed by complete 
dominance QTL, and the lowest for additive QTL. When 
alternate alleles were segregating in the grandparents, the 
power to detect QTL in the POD models decreased for all 
non-POD QTL, except for additive QTL, for which there 
was a slight increase in power compared to having 
homogeneous distributions of alternate alleles in the 
parental breeds (Table 3). 

When the series of POD tests (Tests 2 and 3) were 
applied, most of the non-POD QTL with large or medium 
effects (close to 100%) were not classified as POD under 
the assumption of fixed alleles in the grand parents. The 
accuracy of non-POD declaration was the highest for 
Mendelian QTL across all QTL sizes and mating designs, 
but the lowest for small imprinted QTL, e.g. 24% of 
maternally expressed QTL were declared to be POD in 
design I when alternate alleles were fixed in the parental 
breeds. Spurious POD declaration for imprinted QTL was 
higher in design I than in design II, especially for small 
QTL. When the number of F2 progeny increased, the 
proportion of false POD declaration for imprinted QTL 
decreased, which was more pronounced in design I than in 
design II, when QTL size was not small, or when alternate 
alleles were fixed in the parental Unes. For example, there 
was a 19% and 22% increase of non-POD declaration for 
paternally and maternally expressed small QTL, 
respectively, when parental breeds were fixed for alternate 
alleles in design I (Table 3). 

DISCUSSION 

Detection of polar overdominance QTL 

A least squares regression interval mapping model was 
presented to detect QTL with a unique mode of gene 
expression, POD, in F2 crosses between outbred lines, and a 
series of tests that differentiate POD QTL from Mendehan 
or complete imprinting QTL were developed. Our 
simulation results showed that, in general, POD QTL were 
well detected and characterized when QTL alleles were 
fixed for alternate alleles in the parental breeds (more than 
90% (80%) for medium (small) QTL in F2 mating designs 
that have been used in pig populations). However, when the 
QTL effect was small and QTL alleles segregated in the 
parental Unes, the proposed model and the series of tests did 
not enable effective POD QTL detection, i.e. POD 
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Table 3. Power to detect Mendelian or parent-of-origin QTL under different genetic models, and proportion of the non- polar 
overdominance (POD) QTL declared as POD type based on simulated data 



QTL effect 
(aUele 
frequency)" 



Power to detect QTL 

(%)" 



Declared POD type 



Mend 



Pat Mat 
Expr Expr 



Pod Pod 



QTL effect 
{allele 



Power to detect QTL 

(%)' 



I II Both None frequency)" Mend 



Pat Mat 

Exp Exp 



Pod 
1 



Pod 
II 



Declared POD type 

(%)° 

I II Both None 



Design I (N = 512) 
Pat 0.57 (1.0) 
Pat 0.57 (0.8) 
Pat 0.35 (1.0) 
Pat 0.35 (0.8) 
Pat 0.23 (1.0) 
Pat 0.23 (0.8) 
Mat 0.57 (1.0) 
Mat 0.57 (0.8) 
Mat 0.35 (1.0) 
Mat 0.35 (0.8) 
Mat 0.23 (1.0) 
Mat 0.23 (0.8) 

Design I (N = 1,02^ 
Pat 0.35 (1.0) 
Pat 0.35 (0.8) 
Pat 0.23 (1.0) 
Pat 0.23 (0.8) 
Mat 0.35 (1.0) 
Mat 0.35 (0.8) 
Mat 0.23 (1.0) 
Mat 0.23 (0.8) 

Design II (N = 513) 
Pat 0.57 (1.0) 
Pat 0.57 (0.8) 
Pat 0.35 (1.0) 
Pat 0.35 (0.8) 
Pat 0.23 (1.0) 
Pat 0.23 (0.8) 
Mat 0.57 (1.0) 
Mat 0.57 (0.8) 
Mat 0.35 (1.0) 
Mat 0.35 (0.8) 
Mat 0.23 (1.0) 
Mat 0.23 (0.8) 

Design n (N = 1, 
Pat 0.35 (1.0) 
Pat 0.35 (0.8) 
Pat 0.23 (1.0) 
Pat 0.23 (0.8) 
Mat 0.35 (1.0) 
Mat 0.35 (0.8) 
Mat 0.23 (1.0) 
Mat 0.23 (0.8) 
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Mendelian or parent-of-origin QTL were simulated based on different allele frequencies, magnitudes, and number of F2 progeny in two mating designs I and 11 (see Table 2). 
Add: QTL with additive effect only (d =0), Dom: QTL with complete dominance (a = d), Pat: QTL with paternal expression, JVIat: QTL with maternal expression. QTL effects 
for each type of inheritance were defined, in magnitude order, as large, medium or small QTL, such that the QTL explained 32%, 12.5% or 5.1%, respectively, of the 
phenotypic variance. Alternate QTL alleles were homogenously ( 1 .0/0.0) or differenUy (0.8/0.2) distributed in Fq parental breeds. 

Proportion of rephcates in which a type of QTL was detected at a 5% chromosome-wise level in Mendelian, paternal expression, maternal expression, POD I, or POD II 
models. 

If a type of QTL was detected at a 5% chromosome-wise level in POD models (Testl), a series of POD tests (Tests 2 and 3) were performed for the QTL to be declared as 
POD L POD II, or both types. 



classification down to about 35% for small QTL. When a 
larger number of F2 animals were used, the proportion of 
QTL declared to be POD increased substantially for small 
QTL, e.g. up to 25% for type I POD QTL with segregating 



alleles in the parental lines in design II, where a larger 
number of parents were used. However, there was limited 
gain in the proportion of QTL declared to be POD in design 
I, in which a very small number of parents were used, e.g. 
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only an increase of up to 10% (from 42% to 52%) for type 
II POD QTL with segregating alleles in the parental breeds 
(Table 2). This result could be partly explained by varying 
QTL sizes in realized rephcates of design I. Because few 
animals were used as parents (Fi) and grand-parents (Fq) in 
design I, the realized breed-origin allele frequency 
differences in Fi and F2 generations would be variable 
among the replicates, e.g. more homozygous Fi parents than 
expected, even if they were generated from the parameter 
set of the given allele frequency difference (0.8/0.2), as was 
observed in our previous simulation study (Kim et al., 
2005). This could explain the standard deviation of the 
estimates for POD QTL effects or of the proportion of 
phenotype variance due to POD QTL. Estimates of the 
standard deviation of the QTL variance were larger in 
design I than in design II for POD QTL with the same effect 
and with alleles segregating in the parental lines (Table 2). 
This may prevent some portion of replicates with smaller 
estimates of QTL effect in design I to get through the 
subsequent POD tests (2 and 3) for POD declaration. Also, 
note that proportions of rephcates with POD QTL detection 
in their respective POD models (Test 1) were lower in 
design I than in design II for QTL with segregating alleles 
in the parental breeds (Table 2). 

Our simulation results strongly suggest that POD QTL 
be tested with the presented tests (2 and 3) to differentiate 
POD expression from parent-of-origin expression (complete 
imprinting). When the series of tests of classifying POD 
QTL were not appUed, POD QTL were classified as 
partially expressed QTL, i.e. 100% declaration for large 
QTL with fixation of alternate alleles in the grand parental 
breeds (Table 2). However, when the QTL effect became 
smaller, the proportion of POD QTL declared as partially 
expressed QTL decreased and more POD QTL were 
declared as spurious paternal or maternal expression QTL, 
e.g. up to 35% for small QTL. The rate of the false 
declaration as paternal or maternal expression could be 
reduced by using a larger number of F2 offspring. However, 
the effect of decreasing spurious declaration of complete 
imprinting was in the reverse direction for QTL with 
segregating alleles in the grand parental breeds (Table 2). 

Detection of spurious polar overdominance QTL 

Our simulation results show that most Mendelian QTL 
or complete imprinting QTL with large or medium effect 
were not classified as POD type QTL. When the non-POD 
QTL were analyzed using POD models (Test 1), most of the 
large or medium non-POD QTL were detected in the POD 
models. However, subsequent POD tests (Tests 2 and 3) did 
not allow the QTL to be declared as POD type QTL (Table 
3). However, when the QTL effect was small and breed- 
origin alleles were segregating in the grand parental hues, 
some portion of the imprinted QTL were spuriously 



classified as POD type QTL, which was more pronounced 
in the design (I) with a small number of parents (e.g. Fj 
sires) (Table 3). A similar concern on the use of a small 
number of parents was described in de Koning et al. (2002), 
where spurious detection of imprinted QTL for MendeUan 
QTL was a problem when using an extremely small number 
of Fi sires. Care must be taken, therefore, in the 
interpretation of detected POD QTL with small effects 
when implemented in such a mating design, which includes 
the majority of porcine QTL studies (Bidanel and 
Rothischild, 2002). To decrease false POD classification for 
non POD QTL, use of a larger number of F2 offspring is a 
good option, as shown in our simulation results (Table 3), or 
alternative strategies, i.e. mating designs, or advanced 
lesling methods may be needed. 

It must be noted that detection of imprinting effects, 
including polar overdominance, requires segregation of 
marker alleles within the lines that are crossed, and is not 
possible in inbred Une crosses. Some imprinted genes in 
murines and other mammals, as well as those responsible 
for the callipyge phenotypes in sheep, are clustered in 
chromosomal domains (Tycko and Morison, 2002). Cis and 
trans interactions among imprinted genes around the 
callipyge locus resulted in different meat composition and 
amounts of skeletal muscle between lambs with and without 
callipyge genotypes (Georges et al., 2003). Thus, 
development of advanced imprinting models taking into 
account epistatic relationships between closely linked genes, 
as well as generation of high density maps, may be needed 
to better identify and characterize POD QTL. 

Recently, next generation sequencing technologies have 
enabled the development of high throughput genotyping 
through high density SNP chips such as the 64k Illumina 
Infinium assay, which has enabled genome- wide association 
(GWA) studies to detect QTL with high mapping resolution 
(Lee et al., 2012). Application of a GWA study to detect 
POD QTL may be feasible if the SNP genotypes of Fi 
parents and F2 individuals are known, in which case POD 
tests are based on each SNP with alternate SNP alleles, e.g. 
allele 1 and 2, without the need for genotype information on 
markers that flank QTL, as in this study. 
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